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(54) Post-processing method lor reducing artifacts in block-coded digital Images, and post- 
processing device for actuating such method 



(57) A post-processing method for reducing arti- 
facts in block-coded cfigita) images, characterized by 
providing for: dividing an input image (In) into a plurality 
of image blocks (IB): for each image block (IB), estimat- 
ing global features (GF) of said image block providing 
information on an average content of image edges 
along the horizontal and vertical directions of said 
image block (IB); for each pixel (Px) of an image Mock 
(IB) under examination, estimating local features (LF) 



for said pixel (Px) presiding information on the content 
of image edges along the horizontal and vertical ejec- 
tions of an image area near said pixel (Px); modifying 
the value of said pixel (Px) according to both said global 
features (GF) of the image block (IB) to which said pixel 
(Px) belongs and said local features (LF) of the Image 
area near said pixel (Px). 
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Description 

The present invention relates digital image coding. More precisely, the invention relates to a post-processing 
method for reducing artHacts in Mock-coded digital images, and to a post-processing device suitable for actuating such 
a method. 

With the cfiffusion of digital communication systems, digital images are more and more used. This has led to the 
diffusion of stli and video cameras wilh digital acquisition and processing capability. 

In order to belter exploit storage means and transmission bandwidth, digital image compression standards have 
been developed, such as JPEQ for still images, and MPEG-1 and MPEG-2 for digital television image sequences. 

The above-referred compression standards provides for btock<wfing based on Discrete Cosan e Transform (OCT). 
A digital image is divided into blocks of pixels, and each btock is encoded independently from the others. DCT coeffi- 
cients for the pixels of each block are evaluated and a quantization matrix is applied to the OCT coefficients to reduce 
the information to be stored or transmitted. When the image Is to be displayed, it must be decoded in advance. 

Due to the quantization process, these image compression methods are lossy, I.& they cause a toss of information 
in the decoded image with respect to the original image. The decoded Image can thus present noticeable degradation, 
mainly consisting of two kinds of artifacts known in the art under the names of "grid noise* and "staircase noise". 

In order to reduce the image degradation, post-processing methods of the decoded image have been proposed 
which allow for attenuating grid noise and staircase noise. 

In view of the state of the art described, it ts an object of the present invention to provide a new post-processing 
method for reducing artifacts in block-cod ed tfgitaJ images. 

According to the present Invention, such object is attained by means of a post-processing method for reducing arti- 
facts in block-coded digital images, characterized by providing for: 

a) dividing en input image into a plurality of image blocks; 

b) for each image block, estimati ng global features of said image btock provkfi ng information on an average content 
of image edges along the horizontal and vertical directions of said image block; 

c) for each pixel of an Image block under examination, estimating local features for said pixel providing information 
on the content of Image edges along the horizontal and vertical directions of an image area around said pixel: 

d) modffying the value of said pixel accenting to both said gfobal features of the image block to which said pixel 
belongs and sad local features of the image area around said pixel. 

Aiso according to the invention, there rs aJso provided a post-processing device for actuating said method, charac- 
terized by conprising: 

• first means supplied with an input image for esti ma ting global features of an image block under examination, said 
global features providing information on an average content of image edges along the horizontal and vertical direc- 
tions of said image block; 

- second means supplied with said input image for estimating local features for each pixel of the image block under 
examination, said local features proving information on the content of image edges along the horizontal and ver- 
tical directions of an image area around said pixel; 

• third means supplied with said global featores ami said locsJ features form 
to both said global features arid said local features. 

The features and advantages erf the present invention will be made apparent from the following detailed description 
of an embodiment thereof, illustrated as a non-limiting example in the annexed drawings* wherein: 

Figure 1 is a schematic block diagram illustrating the principle of operation of a method according to the present 
invention; 

Figure 2 shows a digital image divided into image blocks; 

Figure 3 shows in detafl an image bloc* into which the digital image of Figure 2 is divided; 
Figure 4 shows an array of pixels of the image block of Figure 3 ; 

Ftgure 5 shows an image sub-Hoc* of the image btock of Figure 3 used for evaluating global features of the imaQe 
block; 

Figure 6 shows an horizontal processing windewv used tor evaluating focal features in the horizontal direction for a 
generic pixel of the image block; 

Figure 7 shows a vertical processing window used for evaluating local features in the vertical direction for said 
generic pixel; 

Figures 8 and 9 shows two membership functions used to perform a fuzzy 
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Figure 10 is a block diagram of a device according to the present invention; 
Figure 1 1 shows the structure of two blocks of the device of Figure 10; and 
Figure 12 is a block diagram of other two Uocte of the device of Figure 10. 

With reference to Figure 1. there is shown a block diagram illustrating the principle of operation of the post-process- 
ing method according to the present invention. An input decoded compressed digital image in is supplied to a Feature 
Extraction block FE. Sock FE provides for analyzing the image to evaluate global and local features thereof. The global 
and local features, respectively OF and LF. of the image In evaluated by block FE are suppfied as inputs to a Fuzzy 
Process block FUZZY which, according to fuzzy rules, determines parameters FA suitable for determining the Kind of 
filtering to be performed, in accordance to the global and local features GF and FL of the image. The parameters FA 
calculated by block FUZZY are supplied to a Filter Composition block FC which according to said parameters FA deter- 
mines the type of filtering to performed out of a set of predefined filters (block FS). Filter parameters FP determined by 
block FC are then supplied to a Processing block PBOC, also supplied directly with the input image In, which performs 
the filtering of the input image In according to the fitter parameters FP to provide a post-processed output image Out. 

If appears that the Wnd of f Otering to be performed on the decoded input image In Is chosen after an estimation of 
the global art local features of the dec 

filtering is performed, to reduce both staircase noise and grid noise. For areas containing fine details (image edges and 
texture), no filtering is performed. Thus, the method according to the present invention provides for performing a non- 
linear adaptive filtering on the pixels of the decoded image 

The principle of operation outlined above will be now descrfred in dotal. 

As shown in Figure 2, the input image In is partitioned into knage Uocte IB, each containing an equal number of 
pixels. A typical dimension of the blocks Is 8*8 pixels (Figure 3), but this is not however to be intended as a limitation, 
since other block dimensions are suitable 

The image blocks IB of the input image In are scanned Urie by 6ne starting from tte 
one. For each image block IB, the Feature Extraction block FE in Figure 1 determines the global and focal features GF 
and IF. 

Global features of the Image block IB under examination are determined by applying horizontal and vertical Sobe! 
operators: 



to each pixel belonging to an image sub-fclock internal to the Image block IB. For example, the following Sobe! opera- 
torse 



are applied to a 6*6 pixel image sub-btock ISB (gray area in Figure 5). As shown in Figure 4, for each pixel Pi of the 
image suM*ock iSB a 3*3 array of naghbomg pixels M centered in pixel Pi is considered, and the values of the pixels 
of said array M are multiplied by the coefficients of the horizontal and vertical Sobel operators, to obtain: 



horizontal Sobel operator (Hsob): I h21 h22 h33 

Lh31 h32h33. 




rv11 v12 v13i 

vertical Sobel operator (Vsob) :|v21 v22 v23| 

Lv31 v32 v33J 





Hsob - (P3+P8+TP5) • (P1+P6+2T4), 
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Vsob s (P6+P*tf # P7) - (P1+P3t2*P2). 

where P1-P4 and P5-P8 are the values of the pixels (gray levels). 

As known, horizontal and vertical Sobel operators perform a filtering capable of detecting edges along the horizon* 
tal and vertical direction, respectively. 

The output values of the horizontal Sobel operators calculated tor the pixels of image sub-tatock ISB are aocumu- 
laled to obtain an accumulated value Aoc(Hsob), and the output values of the vertical SobeJ operators calculated for the 
pixels of image sub-block ISB are accumulated to obtain an accumulated value Acc(Vsob). Acc(Hsob) gives the high- 
frequency content in the horizontal direction (vertical edges) of the image block IB. Aoc(Vsob) gives the high-frequency 
content in the vertical direction (horizontal edges) of image Hoc* IB. Thus. Acc(Hsob) and Acc(Vsob) respectively pro- 
vide the degree of 'edgeness' of the image biock under examination in the vertical and horizontal directions, rt is to be 
noted that in order to evaluate the global features GF of the image biock IB under examination, only the pixels belonging 
to this block are oonskieied (by applying 3^ Sobel operators to the 6*6 image sut^block ISB, it is not necessary to take 
into consideration pixels belonging to neighboring image blocks). 

Global features GF erf the irriage block under e 
Acc(Vsob). Alternatively, the global features GF of the image block can be formed by an average of the accumulated 
values Acc(Hsob) and Acc(Vscb), to obtain the averag e number of edges in the horizontal and vertical directions in the 
image block under examination. 

Local features LF of the image block IB are estimated according to the following method. All the pixels of the image 
block IB under examination are scanned Dne by line starting from the top-left pixel down to the bottom-right one. To each 
pixel a horizontal processing window is applied: a presented number of pixels respectively preceding and following the 
pixel under consideration and belonging to the same image fine of the pixel under consideration is considered. A surta- 
bl e horizontal processing window HP W is shown in Figure 6. which is a horizontal 1 *5 processing window: tor a given 
pixel, the two preceding pixels Pa. Pb and the two following pixels Pc, Pd belonging to the same Bne are considered In 
Figure 6 there is shown by way of example the horizontal prrxessmg window HPW a^ 
image block. It should be noted that nrtrjn^ the pixels of the irnageblodt IB under ex 

pixels belonging to neighboring image blocks: this is for example the case of the first second, seventh and eight pixel 
of each tine of pixels of the image block IB under examination. 

The horizontal Sob el operator Hsob previously mentianed is applied to each pixel Pa.Pb.Px.Pc.Pdin the horizon- 
tal processing winder HPW, to obtain five output values HS 144S5. Values HS1-HS5 provide the local features in the 
horizontal direction for the pixel under examination Px, i.e. the high-frequency content in the horizontal direction of the 
image region around the pixel under examination. 

Similarly, a vertical processing window is applied to each pixel of the image block IB. The vertical processing win- 
dow is formed by the pixel under consideration Pi and a prescrM 

and preceding and following the pixel under consideration; for example, as shown in Figure 7 the vertical processing 
window VPW can have dimensions identical to the horizontal processing window HPW (5*1 ), and thus contains two pix- 
els Pe. Pf preceding pixel Px and two pixels Pg, Ph following pixel Px in the vertical cfirection. 

The vertical SobeJ operator Vsob premise mentioned sften applied to 
dow VPW to obtain five output values VS1 -VS5. Values VS1 -VS5 form the local features in the vertical direction for the 
pixel under examination, i.e. the high-frequency content in the vertical direction of an image region around the pixel 
under examination. 

The global features GF for the image biock IB under examination (i.e.. the two accumulated values Acc(Hsob) and 
Acc(Vsob) or. in alternative, the average value of Acc(Hsob) and Aoc(Vsob)) and the focal features LF for the pixel 
under examination inside said image block (the ten values HS144S5 and VS1-VS5) are then suppGed to the Fuzzy 
Process block FUZZY. This block provides tor evaluating the degree of membership of a generic value HSi and Vsi 
(U1..5) to two fuzzy sets "Small" and "Big"; these degree of membership can be evaluated by applying to HSi. VSi the 
membership functions depkrted in Figures 8 and 9. In these figures, Thl and Th2 are values depending on the global 
features OF of the image block under examination, Le, on the accumulated values Acc(Hsob) and Acc(Vsob)oronthe 
average of the accumulated values; in the first case, Th1 and Th2 are different for the Hsiand Vsi values; in the second 
case, Thl and Th2 are the same tor Hsi and Vsi values. 

Fuzzy rules having as antecedents the degrees of membership of the output values HSi and VSi to the two fuzzy 
sets "Small' and "Big" are then evaluated. This means that 32 rules are to be evaluated tor both the horizontal and ver- 
tical directions. However, all those fuzzy rules having the same consequence are synthesized in one rule only by an else 
operator. In this way. the system complexity is reduced, and a total of nine rules for each direction have to be evaluated. 

The following fuzzy rules are applied to the five values HS1-HS5 associated to the horizontal oirection: 

1. 0HS1 is Small and HS2 is Small and HS3 e Small and HS4 is Small and HSS is Small, then at ts Big; 

2. //HSI is Small andHS2 is Small andHS& is Small and HS4 is Small and HSS is Big, then o2 b Big; 
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3. H HS1 is Small and HS2 s Small andHS3 is SmaD andHSAis Big andHSS is Small, fhen <x3 csBig; 

4. // HS1 & Small andHSZis SmaO andHSZis Small and HS4s Big and HS5fc Big. Jhen a4fc Big; 
S//HS1 & Small and HS2 is Big andHSQ is Small andHSA is Small andHSSis Small, then aS & Big; 

6. // WS1 is Big and HS2 is Small and WS3 & Small andHS4 is Small andHSSb Small, tf»n a6 & Big; 

7. WS1 cs Big and WS2 is Small and HS3 is SmaJl and HS4 is Small and HS5 is Big, then a7 « Big; 
a tf HS1 is Big and HS2 is Big and HS3 is SmaD and HS4 is Small and HS5 is Small, then aS is Big. 

The activation level of each nd8 depends on the degree of memberships of the pattern of output values HSiof the 
horizontal Sobel operator applied to the five pixels of the horizontal processing window HPW; the degree of member- 
ships depend in turn on the global features GF of the Image block to which the pixel under examination belongs. The 
activation level of the else (ninth) rule is computed as a else « (1 - oave) , where oave is the average activation degree 
of fuzzy rules 1 to 8. at toa8andae!se,arYJasto!arse1of niro 
VS1-VS5) form the output FA of the fuzzy process block FUZZY in Figure 1. 

Each one of the above-listed rules is associated to a respective set of predefined fitter parameters, which are stored 
as a look-up table in block FS of Figure 1 . Suitable predefined filter parameter sets are for example: 

Rule 1: (d 1°1 .0, C12-1 .0, c13°1.0, d 4=1.0. c15=1 .0) if the pixel under examination Oes outside the Image sub- 
block ISB t and (d 1=0.0. c12=1.0, c13«1 A d 4-1 .0, d 5=0.0) if the pixel under examination lies inside the image 
subbiocklSB; 

Rule 2: (c21«05. C22-1.0. C23-1.0, C24-1.0, C25-0.0); 
Rule 3: (031*0.5. C32-1.0. C33-1.0. C34-0.0, c35=0.0); 
Rule 4: (c41=0.5, c42«1.0. C43-1.0. 044-0.0, C45-0.0); 
Rule 5: (c51=0.0, C52-0.0, CS3-1.0. c54=1.0, c55*0.5); 
Rule 6: (c61=0A c62=1.0, C63-1.0, o64=1.0, c65=0.5); 
Rule 7: (c71=0.0. C72.1.0. c73»1.0, C74-1.0, C75-0.0); 
Rule 8: (081=0.0, C82-0.0, C83-1.0, C84-1.0, C85-0.5); 
Else rule: (c91«0.0. C92-0.0. C93-1 X). c94=0.0. c95=0.0). 

The parameters FP of the filter to be applied to the five pixels of the horizontal processing window HPW are calcu- 
lated as a weighted average of the nine fitters described above, with weight factors formed by the activation degrees a1 
to a8 and oteise of the respective fuzzy rules. 

Assuming that oils the actrvat^ 
(ad - aelse ). and oj are the coefficients of the Mh filter (i- 1 ..g, j«i ..5). the weight factor applied to the Mh f flter. asso- 
ciated to the Mh fuzzy rule is: 

R m ai • cq 

and the coefficients Hj of the final horizontal f ater to be applied to the pixels of the horizontal processing window HPW 
are given by: 



£ai-cij 



where N is a normalization factor. 

The horizontally-filtered value Ex of the pixel Px under examination (at the center of the horizontal processing win- 
dow) is then c^latiri as a weigtted averaged the 
zontal processing window HPW, with weight factors formed by the coefficients Hj: 

Px - HfPa* H2*Pb+ H3*Px+ H4'Po H5*Pd. 

Similar calcuJatkxts are performed for the vertical direction, starting from the output values VS1 -VS5 of the vertical 
Sobel operators applied to the pixels Pe, Pf, Px. Pg and Ph in the vertical processing window 
Q-1..5) of the Hter for the vertical direction are calculated in a way orjr^ 
coefficients Hj: 
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9 

VJ N 

5 

where ft ..9) are the activation degrees of nine fuzzy rules for the vertical direction (similar to those feted above for 
the horizontal direction) and cij Q-1..5) now are the predefined filter parameters associated to the Mh fuzzy rule for the 
vertical direction. The coefficients Vj are then applied to the pixels in the vertical processing window VPW to calculate 
a weighted average of the same. The tittered value of the pixel Px under examination, filtered in both the horizontal and 
10 vertical direction, is provided at the output Out of the processing block PROC. 

The value of the pixel Px under examination to be rmdt^ed by the vertcal fitter coefficient 
obtained after having applied to the pixels in the horizontal processing window HPW the horizontal filter Hj f>1 .5): 

Out - VI *Pe ♦ V2*Pf + V3*Px + V4*Pg ♦ VS'Ph. 

iff 

Alternatively, it is possible to evaluate first the vertically-filtered value of the pixel under examination: 

Px « V1 *Pe ♦ V2*Pf + V3*Px ♦ V4*Pg ♦ V5*Ph. 

20 and then performing the fitering in the horizontal direction applying to this value the respective coefficient H3 of the hor- 
izontal filter Hj: 

Out = H1 *Pa + H2*Pb + H3*Px + H4*Pc + H5*Pd. 

25- The sequence is of no importance, the important thing to be underlined being that at the end of the process the value 
of the pixel under examination is the result of both an horizontal and a vertical fatering. 

Figure 10 is a block diagram of a device suitable for actuating the method previously described. The device com- 
prises two main blocks: a first block 1 evaluates the global features QF of the image blocks IB the image to be post- 
processed is divided in, and a second block 2 evaluates the local features LF of the pixels of the image and performs 
30 the filtering acconcSng to both the global features and the local features. 

It is assumed that the image to be post-processed is scanned line by line in a sequential order. Signal In is a stream 
of pixels of the input image scanned line by line. Block 1 is supplied with signal In; signal In also supplies a cascade of 
two One memories LM1 and LM2 whose outputs supplies block 1. 
Inside block 1. signal In and the oulputscif line memories LM1 art 
as implementing a 3*3 pixel window which is used to calculate horizontal and vertical Sobel operators for the pixels of the 
6*6 image sub-block ISB inside each image block IB. Blocks supplies a block 4 which calculates the outputs Hsoband 
Vsob of the horizontal and vertical Sobel operators for those pix^ of the current tr^ 

sub-Wocks ISB of each image block IB. The outputs Hsob and Vsob of block 4 are supplied to an accumulator block 5 
wherein they are accumulated. After eight image Ones. i.e. a fine of image blocks IB. have been scanned, the accumu- 
40 lated values Acc(Hsob). Acc(Vsob) (or alternatively the average thereof) for each image block IB are stored in a mem- 
ory block 6. 

The output of line memory LM2 supplies a cascade of eight further line memories LM3-LM10. Block 2 is supplied 
in parallel with the outputs of line memories LM4-LM10. In this way, evaluation of the local features and calculation of 
the filter parameters starts after block 1 has estimated the global features QF for a tine of image blocks IB. 
45 Inside block 2, a block 7 of pixel delays is suppOed with the outputs of line memories LM4-LM10; by means of the 
line memories LM4-LM10 and the pixel delay block 7 it is poss&fle to implement the 5*1 vertical processing window 
VPW. The outputs L4-L1 0 of the pixel delays block 7 supply a block 8 which applies the vertical Sobel operator to each 
pixel inside the vertical processing window V 

ferred provkfing for calculating five vertical Sobel operators in parallel; the outputs of the five vertical Sobel operators 
50 VS1 -VS5 are supplied to a vertical fuzzy filter 9. which is also supplied with the outputs L6-L10 of the pixel delay block 
7 and the output MOUT of the memory MEM of block 1. MOUT supplies the global features GF of the image block IB 
currently processed by block 2. i.e. the accumulated value Aoc(Vsob) or. alternatively, the average of Acc{Vsob) and 
Acc<Hsob). The vertical fuzzy filter block 9 evaluates the degree of membership of values VS1 -VS5 to the fuzzy sets 
•Small* and "Big" taking into account the global features provided by MOUT, evaluates the activation levels of the nine 
55 fuzzy rules for the verticaJ direction, c^cuJatesthecoeffidents^ 0 D 1..5) of the verti^ f^ and appfies the vertical filter 
coefficients Vj to the five pixels Pe, Pf, Px, Pg. Ph in the verticaJ processing window VPW, to calculate the vertically- 
filtered value gx of the pixel In the middle of the vertical processing window. The output of the vertical fuzzy rater block 
9 forms the verticafly-fStered value Ed of pixel Pd In the horizontal processing window HPW shown in Figure 6 and sup- 
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pGes directly a horizontal fuzzy filter block 10; the output Eg of block 9 also supplies a cascade of four pixel delays D 
whose outputs respectively form the vertically-littered values E& gfe Ea of tfte pixels Pc; Px, Pb, Pa in the horizontal 
processing window HPW and apply the horizontal fuzzy fBter block 10. 

In parallel to blocks 8 and 9. the outputs L7-L9 of the pixel delay block 7 supply a block 11 which applies the hori- 
zontal Sobe! operators to the pixels inside the horizontal processing window HPW. Differently from the vertical sobel 
operators, only one horizontal sobel operator is calculated at a time; a condensation delay block 12 introduces a delay 
for compensating the processing delay of the vertical fuzzy filter block 9. The output of block 12. forming the output of 
the horizontal Sobel operator HS5 applied to pixel Pd of the horizontal processing window in Rgure 6, supplies the hor- 
izontal fuzzy fitter block 10 and a cascade of four pixel delays D, the outputs thereof forming the values HS4, HS3.HS2 
andHSI arri suppling the horizontal fus^ The horizontal fuzzy filter block 10. which is also supplied by 

the output MOUT of the memory block MEM in block 1 providing the value Acc(Hsob) (or alternatively the average of 
values Acc(Hsob) and Acc(Vsob)). evaluates the degree of membership of values HS 1 -HS5 to the fuzzy sets "Small* 
and "Big" according to the value of the global features GF provided by MOUT, evaluates the activation levels of the nine 
fuzzy rules described above tor the filtering in the horizontal direction, calculates the coefficients Hj of the horizontal f fl- 
ter and applies tte parameters Hj to Pb. Px. Pc, Pd in 

the rwrizorttal processing window HPW to obtain the horizontal*- and verticaly-fWered value Out of the pixel Px under 
examination. 

A control circuit CTRL controls the operation of blocks 1 . 2 and the One memories LM1-LM10. 

Rgure 1 1 shows the structure of the vertical and horizontal Sobel operator blocks 8 and 1 1 of Rgure 1 0. They are 
composed in a straightforward way by adders. 

Rgure 12 shows the structure of both the vertical fuzzy fitter block 9 and the horizon 
are the vertical or, respectively, horizontal Sobel operator outputs VS1-VS5 and HS1-HS5. X1-X5 are supplied to a 
fuzzy rule evaluation block 13 which evaluates the activation degrees pi -09 of the nine fuzzy rules for the vertical direc- 
tion or. respectively, the activation degrees a1-a9 of the nine fuzzy rules for the horizontal direction. The activation 
degrees evaluated by block 13 are supplied to a look-up table of respective predefined filter parameters F1 -F9 (forming 
block FS in Figure 1). and the outputs of the look-up table. La the predefined fater parameters crj multiplied by the acti- 
vation degree of the respective fuzzy rule, are supplied to a fater composition block 1 4 which calculates the coefficients 
V1-V5 or. respectively. H1-H5, of the vertical or. respectively, horizontal filter. Said coefficients are then supplied to a 
processing block 15 which is also supplW 

ure 10); block 1 5 applies the fitter coefficients to the pixel values to obtain the filtered value of the pixel under examina- 
tion Px 

Claims 

1. Post-processing method for reducing artifacts in block-coded digital images, characterized by providing tor: 

a) dividing an input image (In) Into a pluraffiy of image blocks (IB); 

b) tor each image block (IB), estimating global features (GF) of said image block providing information on an 
average content of image edges along the horizontal and vertical directions of said image block (IB); 

c) for each pixel (Px) of an image block (IB) under examination, estimating local features (LF) tor said pixel (Px) 
providing information on the content of image edges along the horizontal and vertical directions of an image 
area near said pixel (Px); 

d) modifying the value of said pixel (Px) according to both said global features (GF) of the image block (IB) to 
which said pixel (Px) belongs and said local features (LF) of the image area near said pixel (Px). 

2. Method according to claim 1. characterized in that step d) provides tor: 

d1) defining a set of predefined local features; 

d2) determining degrees of coincidence (a 1 -a9.p1 -p9)of said local features (LF) of the image area around said 
pixel (Px) with each predefined local features of said set. said degrees of coincidence (a1-o9.p1 -p9) depend- 
ing on said global features (OF) of the image block (IB) to which said pixel (Px) belongs: 
<J3) maWng the value of said pixel (Px) equal to a weighted average (E*Ex; of the value of said pixel (Px) and 
of the values of neighboring pixels <Pa-Pd,Pe-Pf), with weight factors (H1-H5.V1-V5) depending on said 
degrees of coincidence of said local features (U^ wfth each d said predefined local features. 

a Method aaxxding to claim 2. characterized in that said determining the degrees of coincidence (a1-«z9,pi-p9) in 
stepd2) provides for performing a fuzzy calculation. 
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4. Method according to daim 3, characterized in that each of said predefined local features is associated to a respec- 
tive group of predefined weight factors (cfl), and each of said weight factors (H1-H5,V1-V5) is calculated as a 
weighted average of corresponding predefined weight factors (dj) of said groups with weight coefficients (a1- 
a9.pl -09} depending on said degress of coincidence of said local features (LP) with each of said predefined local 
features. 

& Method according to claim 4. characterized in that said estimating global features (GF) of the image block (18) 
under examination provides for applying horizontal and vertical Sobel operators (Hsob. Vsob) to pixels belonging to 
an image sub-biock (1S6) internal to said image block (IB) under examination. 

6. Method according to claim 5. characterized in that said estimating global features (GF) of the image block (IB) 
under examination provides for adding outputs of the horizontal Sobel operators (Hsob) applied to each pixel of 
said image sub-block (ISB) to obtain an accumulated output of horizontal Sobel operators (Acc(Hsob)), and adding 
outputs of the vertical Sobel operators (Vsob) appiied to each pixel of sati 

mulated output of vertical Sobel operators (AccfVsob)). 

7. Method according to claim 6. characterized in that said global features (OF) of the image block (IB) under exami- 
nation are formed by said accumulated outputs of the horizontal and vertical Sobel operators 
(Acc{Hsob)AxHVsob)). 

8. Method according to claim 6, characterized in that said global features (GF) of the image block (IB) under exami- 
nation are formed by an average of said accumulated outputs of the horizontal and vertical Sobel operators 
(Acc(Hsob)AxXVsob)). 

9. Method according to claim 7 or 8, characterized in that said estimating focal features (LF) for a pixel (Px) of the 
image block (IB) under examination provides for: 

d) considering an horizontal processing window (HPW) containing the pixel (Px)' under examination and 
neighboring pixels (Pa.Pb.Pc.Pd) belonging to a same image Brie as the pixel (Px) and preceding and following 
the pixel (Px): 

c2) applying said horizontal Sobel operator (Hsob) to each pixel (Pa-Pd.Px) of the horizontal processing win- 
dow (HPW) to obtain a horizontal pattern of horizontal Sobel operator outputs (HS1-HS5); 
c3) considering a vertical processing window (VPW) containing the pixel (Px) under exarrtinatipn and neighbor- 
ing pixels (Pe.Pf.Pg, Ph) belonging to a same column of pixels as the pixel (Px) and preceding and following the 
pixel (Px): 

c4) applying said vertical Sobel operator (Vsob) to each pixel (Pe-Ph,Px) of the vertical processing window 
(VPW) to obtain a vertical pattern of vertical Sobel operator outputs (VS1 -VS5). 

10. Method according to claim 9. characterized in that said horizontal and vertical processing windows (HPW.VPW) 
contains each orie f ive pixels and are centered at said pixel (Px) under examination. 

11. Method according to claim 10. characterized in that step d2) provides for determining degrees of membership of 
each rtorizontaJ Sobel operator output (HS1 -HS5) of the horizontal pattern to a first fuzzy set "Small" and to a first 
fuzzy set "Kg", evaluating activation degrees (al -a9) of a f irst set of fuzzy rules each one associated with at least 
one predefined horizontal pattern of horizontal Sobel operator outputs, determining degrees of membership of 
each vertical Sobel operator output (VS1-VS5) of the vertical pattern to a second fuzzy set "Small" and to a second 
fuzzy set "Big' and evaluating activation degrees (p1-£9) of a second set of fuzzy rules each one associated with 
at least one predefined vertical pattern of vertical Sc^ operator outputs, 

12. Method according to claim 11, characterized in that said determining degrees of membership of the horizontal 
Sobel operator outputs (HS1-HS5) of the horizontal pattern to said first fuzzy sets "Small" and "Big" provides tor 
determining a first and a second membership functions depending on said global features (GF), and said determin- 
ing degrees of membership of the vertical Sobel operator outputs (VS1-VS5) of the vertical pattern to said second 
fuzzy sets "SmalT and "Big" provxtes for cUrtermming a third and fourth membership functions depending on said 
global features (OF). 

13. Method according to daim 1 2. characterized in that said groups of predefined weight factors (dj) comprise groups 
of predefined hcwizontal weight factors and groups of predefined vertical weight factors, each fuzzy rule of said first 
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86t being assodated with a respective one of said groups of predefined horizontal weight factors, and each fuzzy 
rule of said second set being associated with a respective one of said groups of predefined vertical weight factors. 

14. Method according to daim 13. characterized in that said weight factors (H1 -H5.V1 -V5) comprise horizontal weight 
factors (H1-H5) and vertical weight factors (V1-V5). said horizontal weight factors (H1-H5) being determined by 
making a weighted average of the predefined horizontal weight factors with weight coefficients being formed by the 
activation degrees (a1-a9) of the fuzzy rules of the first set and said vertical weight factors (V1-V5) being deter- 
mined by making a weighted average of the predefined vertical weight factors with weight coefficients formed by 
the activation degrees (01 -p9) of the fuzzy rules of the second set 

1 5. Method according to daim 14, characterized in that the value of the pixel (Px) under examination is rmxfified by 
applying the horizontal weight factors (H1-H5) to the values of the pixels (Pa~Pd.Px) in the horizontal processing 
window (HPW) and applying the vertical weight factors (V1 -V5) to the values of the pixels (Pe-Ph.Px) in the vertical 
processing window (VPW). 

16. Post-processing device for reducing artifacts in block-coded digital images, characterized by comprising: 

- first means (1) supplied with an input image (bi) for estimating global features (MOLTT) of an image block (IB) 
under examination, said global features providing information on an average content of image edges along the 
horizontal and vertical directions of said image block (IB); 

- second means (7,8,1 1 ,12,D) supplied with said input image (In) for estimating focal features (VS1-VS5.HS1- 
HS5) for each pixel (Px) ol the image Hock (IB) under examination, said local features providing information on 
the content of image edges along the horizontal and vertical directions of an image area around said pixel (Px) ; 

• third means (9.0110) supplied with said global features (MOUT) and said local features (VS1-VS5.HS1 -HS5) 
for modifying the value of said pixel (Px) according to both said global features (MOUT) and said local features 
(VS1-VS5.HS1-HS5). 

17. Device according to daim 1 6, characterized in that said first means (1) comprises means (3,4) for evaluating hori- 
zontal and vertical Sobel operator outputs {Hsob, Vsob) of horizontal and vertical Sobel operators appfied to pixels 
of an image sub-block (ISB) internal to said image block (IB) under examination, and accumulator means (5) for 
accumulating the horizontal Sobel operator outputs (Hsob) and the vertical Sobel operator outputs (Vsob) for each 
pixel of the image sub-block (ISB). 

18. Device according to daim 17, characterized in that satf second 

for evaluating vertical Sobel operator outputs (VS1-VS5) of vertical Sobel operators (Vsob) applied to said pixel 
(Px) and to verticafly^eighboring pixels (Pe.Pf.Pg.Ph) prececfng and following the pixel (Px) in the vertical direc- 
tion, and fifth means (11.12.D) for evaluating horizontal Sobel operator outputs (HS1-HS5) of horizontal Sobel 
operators (Hsob) applied to said pixel (Px) and to horizorrtailvHieighboring pixels (Pa.Pb.Pc.Pd) preceding and fol- 
lowing the pixel (Px) in the horizontal direction. 

1 9. Device according to claim 18. characterized in that said third means (9.D.1 0) comprises vertical filter means (9.D) 
supplied with satf vertical Sobd 

ticalty-fOtered value (M of said pixel (Px) depencfing on said global features (MOUT), said vertical Sobel operator 
outputs (VS1-VS5) and the values of said verttaUy-neigrtboring pixels (Pe.Pf.Pg.Ph), and horizontal filter means 
(10) supplied with said global features (MOUT), said horizontal Sobel operator outputs (HS1-HS5) and said verti- 
caDHBtered value (gxj of said pixel (Px) for caJcufattng a horizontal 

said global features (MOUT), said horizontal Sobel operator outputs (HS1 -HS5) and the verticafly-ffiered values 
(fefoP&gfcEg) of said pixel (Px) and said horizontally-rteighboring pixels (Pa.Pb.Pc.Pd). 

2a Device according to daim 1 9, characterized in that said vertical f Star means (9.D) comprises first fuzzy computation 
means (13) supplied with said vertical Sobel operator outputs 
uating degrees of coincidence 01-09) of me pattern of verttal Sobel operator o^ 

dsfmad patterns of vertical Sobel operator outputs, said degrees of coincidence 01 -09) depencfing on said global 
features (MOUT), a look-up table of predefined vertical fttter coeffidents (F1-F9) each one associated to at least 
one respective predefined pattern of vertical Sobd operator oufou^ 

(14) supplied with said predefined vertical fater coefficients for generating a set of vertical filter coeffidents (V1-V5) 
which are a weighted average of the predefined vertical filter coeffidents with weight coeffidents formed by said 
degrees of coincidence (P1-p9). and a vertical fflter (15) supplied with said vertical tier coeffidents (V1-V5) and 
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the values of the pixel (Px) and the vertjcaily-naghboring pixels (Pe.Pf, Pg.Ph) for providing at an ouput (FOUT) said 
verticaly-filtered value of the pixel (Px). 

21. Device according to claim 20, characterized In that said horizontal filter means (1 0) composes second fuzzy com- 
putation means (13) supplied with said horizontal Sobel operator outputs (HS1-HS5) and said global features 
(MOUT) for evaluating degrees of coincidence (ai-o9) of the pattern of horizontal Sobel operator outputs (HS1- 
HS5) with a set of predefined patterns of horizontal Sobel operator outputs, said degrees of coincidence (a 1 -a9) 
depending on said global features (MOUT), a look-up table of predefined horizontal filter coefficients (Fl -F9) each 
one associated to at least one respective predefined pattern of horizontal Sobel operator outputs, a horizontal filter 
coefficients composition means (14) supplied with said predefined horizontal filter coefficients (F1-F9) for generat- 
ing a set of horizontal filter coefficients (H1 -H5) which are a weighted average of the predefined horizontal filter 
coefficients with weight coefficients formed by said degrees of coincidence (at -a9). and a horizontal f flter ( 1 5) sup- 
pied with said horizontal filter coefficients (H1-H5) and the vertically-f iltered values (Pg.Pb.Px.Pc Pd) of the pixel 
(Px) and the horizomafly^ghboring pixels (Pa.Pb,Pc.Pd) tor rjroviding at an ouput (FOUT) said horizontaDy-fil- 
tered value (Out) of the pixel (Px). 
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